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Abstract 



> . 

The probability distribution functions (PDFs) of momentum flux and 



zonal flow formation in ion-temperature-gradient (ITG) turbulence 
are investigated in two different models. The first is a general five 

o ; 

field model (rij, 4>, T$, T e , v^) where a reductive perturbation method 
is used to derive dynamical equations for drift waves and a zonal 

•i-H . 

^ | flow. The second is a reduced two-field model (</>, Tj) that has an 

exact non- linear solution (bipolar vortex soliton). In both models the 
exponential tails of the zonal flow PDFs are found with the same 
scaling (PDF ~ exp{— czf^zfW but with different coefficients czf- 
The PDFs of momentum flux is however found to be qualitatively 
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different with the scaling (PDF ~ exp{— cmR s }), where s = 2 and 
s = 3/2 in the five and two field models, respectively. 
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I Introduction 



One of the main challenges in magnetic fusion research has been the pre- 
diction of the turbulent heat and particle transport originating from various 
micro-instabilities. The ion-temperature-gradient (ITG) mode is one of the 
main candidates for causing the anomalous heat transport in core plasmas 
of tokamaks [lj. Significant heat transport can however be mediated by 
coherent structures such as streamers and blobs through the formation of 
avalanche like events of large amplitude, as indicated by recent numerical 
studies [2]- [I] . These events cause the deviation of the probability distribu- 
tion functions (PDFs) from a Gaussian profile on which the traditional mean 
field theory (such as transport coefficients) is based. A crucial question in 
plasma confinement is thus the prediction of the PDFs of the transport due 
to these structures and of their formation [5]- [9]. Note that a renormalized 
perturbative theory can easily make an considerably large error in predicting 
PDFs where the coherent structures are crucial. 

In previous papers the Hasegawa-Mima model of drift wave turbulence 
was used to show that the PDF tails of global momentum flux and heat 
flux are significantly enhanced over the Gaussian prediction j5]- [6]; this was 
later shown to hold also in ITG mode turbulence [Tj- [9]. Specifically, the 
tails of PDF of global momentum flux and heat flux have been shown to 
be stretched exponential with the form ~ exp{— c(R/ Rq) 3 / 2 } [5]- [9], which 
is broader than a Gaussian. These results provided a novel explanation for 
exponential PDF tails of momentum flux found in recent experiments at 
CSDX at UCSD [10]. Note that considerable transport can be mediated 
by rare events of high amplitude at the PDF tails even if the latter have a 
low amplitude. Since the PDF tails of momentum flux are enhanced over 
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the Gaussian prediction in this case, these events are more likely to mediate 
considerable transport. In fact, it was found that for certain values of param- 
eters, large events are crucial for transport. Specifically, the overall amplitude 
was shown to be larger in ITG than in HM turbulence for reversed modon 
speed (U < 0). The PDF tail of zonal flow in ITG turbulence was also sig- 
nificantly increased compared to that in HM turbulence. Furthermore, zonal 
flows are shown to be generated more likely further from marginal stability, 
which will then regulate ITG turbulence, leading to a self-regulating system. 
Namely, while ITG turbulence is a state with high level of heat flux, it also 
generates stronger zonal flows that inhibit transport. This also suggested 
that stronger zonal flows are generated in ITG turbulence compared with 
ETG turbulence. It was also shown that shear flows can significantly reduce 
the PDF tails of Reynolds stress and zonal flow formation [5]- [9]. 

The purpose of the present paper is to investigate the likelihood of the 
formation of coherent structures by computing the PDF (tails) of zonal flow 
formation and the PDFs of momentum flux [llj by using the two different 
ITG turbulence models. The first model is a widely applicable five field 
model (rij, 0, Tj, T e , t>j||) of the ITG turbulence coupled to an ion vorticity 
equation for the zonal flows. For the computation of the PDFs, a reductive 
perturbation method is used to obtain dynamical equations for drift waves 
and zonal flows. Due to the perturbation method, this model is only valid 
for weak drift wave and zonal flow electric potential. Moreover an exact 
non-linear solution in this model is not available that can be used as an 
ansatz for a coherent structure for the intermittent transport. However, the 
computation of the PDF tails requires only mean values over the coherent 
structures. The second is the two field (0, T«) model which has an exact non- 
linear bipolar vortex soliton (modon) solution that will be utilized in the 
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computing PDFs. Note that the two-field fluid model for the ITG mode [12J- 
[T5] has been successful in reproducing both experimental [16] and non-linear 
gyro- kinetic results [T7] . 

The theoretical technique used here is the so-called instanton method, a 
non-perturbative way of calculating the PDF tails. The PDF tail is first for- 
mally expressed in terms of a path integral by utilizing the Gaussian statistics 
of an external forcing with a short correlation time. An optimum path will 
then be associated with the creation of a coherent structure (among all possi- 
ble paths) and the action is evaluated using the saddle-point method. In the 
second model this coherent structure is identified as the modon. The saddle- 
point solution of the dynamical variable 4>(x, t) of the form t) = F(t)ij)(x) 
with F(t) which is localized in time is called an instanton. Since the instan- 
ton exists during the formation of the coherent structure, the bursty event 
can be associated with the creation of a coherent structure. Note that the 
function ip(x) here represents the spatial form of the coherent structure. 

One of our main results is that the PDF tail of the zonal flow formation 
is PDF ~ exp{— czFj4 > %F\i U = 1,2 to denote the different constants in 
the two models), in agreement with earlier findings [9J. Interestingly the 
predicted PDF tails of momentum flux are Gaussian when the feedback of 
zonal flow on the turbulence is incorporated. This result differs from what 
has theoretically been found earlier [9] whereas it is in agreement with pre- 
dictions from non-linear simulations of turbulence [19]- [20] and gyro-kinetic 
toroidal simulations [21]- [22 j . The reason for this is that previously the 
zonal flow has been treated as passively evolving by the turbulence without 
the feedback loop to turbulence. This feedback gives a cubic non-linearity in 
Eq. (7) describing the ITG fluctuations which changes the characteristics of 
the turbulence. Physically, it is because the feedback of zonal flow regulates 
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turbulence, limiting its growth from turbulence. 

The result from both models support the conclusion that while ITG tur- 
bulence maintains high level of transport, this may be suppressed by shear 
flow. Zonal flows are also shown here to have an enhanced likelihood of the 
generation further from marginal stability which will then regulate the ITG 
turbulence (which is more prominent with increased shear flow) leading to a 
self-regulating system. 

The paper is organized as follows. In Sec. II the drift wave - zonal flow 
system is presented and in Sec. Ill and IV the computation of the PDF and 
a numerical study of the five field model are presented. In Sec. V and VI 
the PDF tails are computed and a numerical study of the reduced two field 
model are outlined. The paper is concluded by Discussion and summary in 
Sec. VII. 

II The drift wave - zonal flow system 

In the five-field model the ITG mode turbulence is modeled using the conti- 
nuity and temperature equations for the ions and assuming the electrons to 
be Boltzmann distributed, closely following previous papers Ref. [12]- |15j . 
The ion continuity, temperature and parallel ion momentum equations then 
become, 





+z/ vV = -[<M + [0,vi0j +f 

dT t ( 2\d<\) 2dn 



(1) 



-[0,T i ] + -[0,n] + 5 Tj , 




(2) 
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Eqs. (l)-(3) are closed by using the quasi-neutrality condition. Here [A, B] = 
(dA/dx)(dB/dy) - (dA/dy) (dB/dx) is the Poisson bracket; v = 0.78 x 
lCT 12 (no/T i o' /2 )(r/.R) is a neoclassical damping; / is a forcing; n = (L n /p s )5n/no, 
<p = (L n /p s )e5(p/T e , Ti = (L n /p s )5Ti/T i0 and v t \\ = (L n /p s )5vi\\/c s are the 
normalized ion particle density, the electrostatic potential, the ion temper- 
ature, and the ion parallel velocity respectively. In equations (1) and (2), 



t = Ti/T e , p s = c s /Q ci (c s = JTe/rrii, Q ci = eB/m,ic). We also used 



L f = - (dlnf/dr) 1 (/ = {n,Tj}), rji = L n /L Ti , e n = 2L n /R where R is 
the major radius and «j = r(l + ?]j). The perpendicular length scale and 
time are normalized by p s and L n /c s , respectively. The geometrical quan- 
tities are calculated in the strong ballooning limit (9 = 0, gi (9 = 0) = 1.0) 
with u± = k y v+ = p s c s k y /L n . In the ion temperature equation [Eq. (2)] we 
have included a heat source Sxi given in Appendix. For the zonal flow Eq. 
(3) vanishes (k\\ = 0) and we cannot use Boltzmann distributed electrons. 
Thus we will employ the following ion vorticity equation [Eq. (4)] and simple 
electron energy equation [Eq. (5)] 



Here r\ e = L n / L Te . The damping (with strength v) represents a neo-classical 
damping of the zonal flow. Note that the damping will affect the full dynam- 
ical system in Eq. (7)- (8). The simple electron energy equation is motivated 
by the fact that we are interested in the modes propagating in the ion drift 
direction. 

In order to find the coupled system of equations for the ITG mode and 
the zonal flow we utilize the reductive perturbation method wherein the 




d d \ ( dn dT 

7^-^--Vijvi0 + ^^(l + .)- + ^ 




= [0,Vi0] (4) 
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perturbed variable a is written 

a(x,y,t) = ££4 +{2/3 ^x{ 1+(2/ ^ 

+ ^ e (4/3) + (2/3)n (T [(4/3) + (2/3)n] (;r ^ )0 
n 

Here the coordinates are scaled as £ = e^ 3 (y — At) and ( = t^ z t with a small 
expansion parameter eo ~ fr ~ 10~ 2 . The drift wave - zonal flow system 
then becomes 

Ca^ + iC 2 ^ + C 3 o 0i = -wC^ + f, (7) 

Here, the coefficients Cj, C 2 , C 3 , C 4 , D 1; £> 2 and _D 3 are complex numbers 
whose forms are provided in the appendix; the variables with tilde and bar 
denote average over x. Note that if the time derivatives are neglected in Eqs. 
(7)- (8), a non-linear Schrodinger (NLS) equation for the fluctuating potential 
is obtained |15|. 



Ill Non-perturbative calculation of zonal flow 
PDF in the five field model 

We calculate the PDF tails of momentum flux and zonal flow formation by 
using the instanton method. To this end, the PDF tail is expressed in terms 
of a path integral by utilizing the Gaussian statistics of the forcing / [18J. 
The PDF of Reynolds stress and zonal flow formation can be defined as 

P(Z) = (S(Zo-Z)) 

dX exp(iAZ) (exp(— i\Z )) 

d\exp(i\Z)I x , (9) 
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where 



h = (exp(-iAZo)). 



(10) 



I\ in (9)-(10) can then be rewritten in the form of a path-integral as 



In the following we assume that the zonal flow potential is stationary -J^ = 
on the time-scale of the fluctuations. The effective action S\ in Eq. (11) can 
be then expressed as, 



Recall £ and ( are the scaled spatial and time variables, respectively. To 
obtain Eq. (12) we have assumed the statistics of the forcing / to be Gaussian 
with a short correlation time modeled by the delta function as 



and (/) = 0. The delta correlation in time was chosen for the simplicity of 
the analysis. In the case of a finite correlation time the non-local integral 
equations in time are needed. 

IV The PDF tails and numerical studies for 
the five field model 

We have expressed the PDF in terms of a path-integral. An approximate 
value of this path-integral can be found for large values of the parameter 




(11) 




(12) 



(/(e,o/(r,o> = *(c-c')«(e-o, 



(13) 



9 



A — > oo, by using a saddle point method. The action in Eq. (12) can be 
expressed using the instanton ansatz 0i(£, C) = ipi (0-^(0 an d 0o(£,C) — 
^(OG'(C) as, 



Sx = -i 

«o 
2 



+ i 



/ rfC (CiFFi + iC 2 k\FF x + C 3 GFF 2 + ivC A FF^j 

i\% J d(G(()5(() 
- % J d(F 3 (D 2 G- RD 3 F 2 ). (14) 

Here k is the strength of the forcing k. The parameter k^(E C) is the inverse 
length scale in the £ direction. In Eq. (14) we have defined the variables, 

$ = Jd^o, (15) 

Fi = JdZMu (16) 

/ d£Ml (17) 

/ dfr^o, (18) 

F A = J d£4wl>l = RF z . (19) 



F 2 

^3 



Note that the relation between the conjugate variables F 3 and F 4 is based 
on the fact that the zonal flow is driven by the Reynolds stress and also that 
the parameter R(e C) has to be determined in such a way to make Reynolds 
stress real. Since the saddle point action is determined by the extremum of 
the action, we require the first functional derivatives of the action to vanish; 

^ = -t(-C 1 F l -tC 2 klF 1 + C 3 F 2 G + tvC 4 F 1 ) 

- 2iRD 3 FF 3 = 0, (20) 
6S * _ ;tn t? in 1.2; 



-%{C X F - iC 2 k\F + ivC 4 F) + kqF\ = 0, (21) 

= —iC 3 FG + k F 2 = 0, (22) 

= D 2 G-RD 3 F 2 = 0, (23) 
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SFt 

SSx 
5F 2 

SF 3 



^ = -iC 3 F 2 F-iF 3 D 2 + i\Q 6(<;)=0. (24) 

The initial conditions F(0) = F and G(0) = G are found from Eqs. (23)- 
(24) as, 

- ( ^^xm (25) 



V C§RD 3t 
RD 3 ( iD 2 <S> 



1/2 



The instanton solution is found from Eqs. (20) - (24) for ( < 0. Starting 
with Eq. (20) and then by substituting the conjugate variables in Eqs. (21) 
- (24) we obtain 

F = ^F 2 = Vl F + 3 V2 F 5 , (27) 



where 



Vi = ^ 2 {-C 2 2 k\ + 2C 2 C A vk\-v 2 C 2 ), (28) 



\_ (cjifn 

Cl\ Dl 



V2 = ^ V ' ( 29 ) 



The integration of Eq. (27) in time, gives us the instanton solution as, 



//(C) = H exp{2^(}, (31) 
#o = F * = , (32) 

m + V^l V ^2^0 + »7l 



where F is given in Eq. (25). The next task is to compute the action 
integrals and estimate the A dependence of S\ in the large A — > oo limit. We 
start by substituting Eqs. (20) - (23) into the action in Eq. (14) keeping 
only the highest powers of A, 

,R 2 D 2 „ R \ 



h » i j^dt\ClF 2 + Ci 



3 n 2 

'-oo \ 
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RD 3 n , , 

+ zA$ -^F 2 (33) 



J Q F °dF(2ClF) +i \^FZ 



RD 3 ( iD 2 <$> 



2 
1/2 



ih ZF \ 3/2 . (35) 



Here, 



^ = $0 ^7(-cfi^J • (36) 

The PDF is then found from Eq. (9) by utilizing the saddle point method 
P(Z) = f d\e~ iXZ - s * (37) 



We find a Aq such that /(A) = —iXZ — ihzF^ 3 ^ 2 attains its maximum and 



compute that value. This gives Aq = ( ) and 



4 .3 



2 



/(A »» = ~2Kf F Z ' (39) 

with the PDF 

P{Z) ~ e" ez3 , (40) 
4 

27^ 



© = ^V7T2 • (41) 

The PDF tails found in the five field model [Eq. (41)] has the same expo- 
nential form as in Ref. [9]. Although no exact non-linear solution was used 
as an ansatz for the coherent structure the system, we still obtain the PDFs 
with the same exponential form. To elucidate on the quantitative differences 
the results between the two and five field models we show the PDF tails by 
using different values of temperature gradient (rji) in Figures 1 and 3. 

In Figure 1 the PDF tails of zonal flow formation obtained in the five-field 
model [Eq.(40)] is presented for r\i = 2.0 (red, dashed line), r\i = 4.0 (blue, 
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iog(* ZF ) 



Figure 1: The PDF tail of zonal flow formation in the five field model for 
r/i = 2.0 (red, dashed line), rji = 4.0 (blue, solid line) and rji = 6.0 (black, 
dashed-dotted line). The parameters are r = 2.0, e n = 1.0, g^ = 1, k x = 0.3, 
k y = 0.3, r] e = 0.0, k = 3000. 

solid line) and rji = 6.0 (black, dashed-dotted line). Other parameter values 
are r = 2.0, e n = 1.0, gi = 1, k x = 0.3, k y = 0.3 and r\ e = 0.0, k = 3000. It 
is interesting to notice that the PDF tails of zonal flow formation in Fig. 1 
are suppressed in comparison with the results found in the reduced two field 
model shown in Figure 3A for similar parameter values. 

Note that the usefulness of the predicted PDF of zonal flow formation 
above is limited since the coefficient G in Eq. (40) is dependent on complex 
coefficients Ci, C2, C2, C3, C4, Di, D 2 and D 3 from the original system of 
Eqs. (7) - (8). As is seen in the appendix due to the complexity of the 
coefficients it is rather difficult to give an explicit region where the analysis 
is useful. 

The analysis can easily be extended to find PDFs of the momentum flux 
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in the system of Eqs. (7) - (8), with the action S\ 

S x = -if d( (C 1 FF 1 + iC 2 k\FF x + C 3 GFF 2 + ivC^FF^ 
+ K ijdC(Fl + F\ 



+ z\R J d(F(0 2 5(() 

- % f d(F 3 (D 2 G-RD 3 F 2 ), (42) 



where i? is the Reynolds stress 

= (43) 



Note that an average over x has already been taken and that the momentum 
flux is non-zero since there is a natural phase shift introduced by the electron 
physics. A similar equation for the instanton can be found, but with different 
values of F and G at t — 0, 



Fo = ^VA, (44) 
Go = ^^A, (45) 

The PDF tail of momentum flux can then be computed by following a similar 
analysis with the result 

P(R) ~ exp{-0 M (^) 2 }, (46) 
©m = (47) 



where 



h M = ~Vv~2 + ^. (48) 



1 ^_ 2k 

Interestingly the exponential forms of the PDFs of momentum flux are qual- 
itative different from what was found in ITG turbulence [7j- [9] and in 
Hasegawa-Mima turbulence [3]- [5] which was the stretched PDF ~ exp{— cR 3 ^ 2 }. 
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Specifically the predicted PDF tails of momentum flux are Gaussian when 
zonal flow feeds back on the turbulence. This result differs from what has 
theoretically been found earlier [9] where the zonal flow is treated kinemati- 
cally. It is however in agreement with predictions from non-linear simulations 
of turbulence [12]- [20] and gyro- kinetic toroidal simulations [21]- [22]- The 
reason for this is the cubic non-linearity in Eq. (7) describing the ITG fluctu- 
ations while previously the zonal flow was treated as passively evolving by the 
turbulence. Physically, it is because the feedback of zonal flow regulates tur- 
bulence, limiting its growth from turbulence. It is interesting that although 
our model for turbulence and zonal flows is expected to have a limited region 
of validity due to the perturbation expansion, it captures the main features 
of what has been found in numerical simulation of plasma turbulence with 
zonal flows. 

Note that non-linear models with cubic non-linearities for turbulent fluc- 
tuations was obtained using perturbation theory and are thus valid only for 
small values of the amplitude (small (pi). The statistical property of the fluc- 
tuations differs radically in systems with cubic and quadratic non-linearities 
as shall be discussed in more detail in Sec. VII. The predictive power of the 
method used here should be able to discriminate between models since the 
statistics of the fluctuations and zonal flow may change qualitatively. We 
note that, in ETG turbulence an evolution equation for the fluctuations with 
a cubic nonlinearity has previously been found [23] . 
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V The PDF tails in the reduced model two 
field model 



Due to the perturbation method, the five field model considered in Sec. III-IV 
is only valid for weak drift wave and zonal flow electric potential. Moreover 
an exact non-linear solution in this model is not available that can be used 
as an ansatz for a coherent structure for the intermittent transport. We will 
thus now consider the reduced two field model where a non-linear solution 
(modon) can be found and no reductive perturbation expansion is needed. 
It is a reduced model where the ITG mode turbulence is modeled using the 
continuity and temperature equation for the ions with adiabatic electrons [5]- 
[9]. In this model the effects of parallel ion motion, magnetic shear, trapped 
particles and finite beta on the ITG modes are neglected since they were 
shown to be not critical in previous work. Unlike the five-field model, the 
generation of a zonal flow is treated passively while the background fluctu- 
ations are affected by an imposed sheared velocity Vq. We note that, the 
effects of mean flow on the zonal flow are weak [26]- [27] while the effects of 
mean flow on the turbulence itself is much more prominent. 

We formally calculate the PDF tails of momentum flux and zonal flow 
structure formation by using the instanton method. The probability distri- 
bution function for Reynolds stress Z\ = R or zonal flow structure formation 
Z2 = 4>zf can be defined as in a generalized Eq. (9) where 



The integrand can then be rewritten in the form of a path-integral as 



Here, the parameter j refers to the two specific cases included in the present 



{exp(-iXjZj)). 



(49) 




(50) 
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study. The angular brackets denote the average over the statistics of the 
forcing /. By using the ansatz % = x4>, the effective action S\ j in Eq. (11) 
can be expressed as, 

Sx 3 = - ? /rf 2 ^0^-(^-a^)Vi0 + \/ o (l-V 2 )0 

+ ^ J d 2 xd 2 x'(f>(x)K(x - x')(f>(x') 

+ i\ 2 J d 2 xdt(j)zF${t) 

+ J d 2 xdt(<j) ZF + R (v x v y }). (51) 
InEq. (51), 



P = 1 + - + -X, (52) 

T T 



v - ^-U 1 ~ U + V ^ (53) 
dip dtp\ 



The PDF tails are found by calculating the value of S\ j at the saddle- 
point in the two cases; the PDF tail of momentum flux by taking into account 
the effect of a shear flow (Ai — > oo, A2 = 0) and the PDF tail of structure 
formation of zonal flow (Ai = 0, A 2 — > 00). The integral in Eq (51) is divided 
in five parts; K\ the ITG integral; K 2 the forcing integral; K 3 the momentum 
flux integral; K4 the zonal flow integral; and finally K$ represents the zonal 
flow evolution integral. The action in the first case can be found as; 

S Xl ~ —ihXl (55) 
h = K 1 + K 2 + K 3 + K± + K 5 , (56) 
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... H 1 H a 



+ 24 ( 1 + 6 -)^(3(7^ip-4(7^ip)J' < 58 > 

X 3 = #0^(0), K 4 = 0, (61) 
2iZoVC 1 

= h^~V (62) 

^ = 4A ~ 2 ' A= (TT^- (63) 

Here C = 2kq\\/((1 + e 2 )7), 7 = ci(l + fc 2 + 2a/ k 3 ) (k is the modon wave 
number), ci = aa/Ji(ka) (a is the modon size), a = Ai — k 2 A 2 , Ax = 
(l-e n -U + V Q )/P, A 2 = ( ai + U- V )/P, 72 = k/3a/2, V = (1 + 6e 2 ) 72 2 
and k is the strength of the forcing k(x —y). To obtain the integrals in Eq. 
(51) we have used the localized modon solution, involving the parameter e 
as done in Refs. [8]- [9]. The PDF tail of the Reynolds stress (R) can now 
be found by performing the integration over Ai in Eq. (51) using the saddle- 
point method in the same fashion as done in Ref. [S]- [9] i.e. recall Eq. (9) 
gives P(R) ~ / d\i exp{— i\\R — S\A- ~ J dXi exp{— i\iR + ihXf}. Now, 



the saddle point integral is evaluated at the saddle point Ximax = yR/(3h) 
which maximizes P(R) with the result 

P(R) ~ exp{-^A 3/2 }, (64) 

* = Iwh (65) 



This result is similar to the previous ones where a similar exponential PDF 
was found [5]- [9]. However, here the focus will be on how the coefficient £1 
depend on the imposed shear flow (Vq). 
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In the second limit where the limit of Ai = and A 2 — > oo gives the PDF 
tail of the structure formation itself, the action can be shown to be 

S X2 ~ -ih\ 3 2 /2 , (66) 

with 

K 3 = 0, (67) 



= 2<j) ZF0 ^-, (68) 
C = 2/ t0 A 2 /((l + e 2 ) 7 ). (69) 

The PDF tail can now be computed by following similar analysis, with the 
result 

P^zf) ~ exp{-6(0zF/0ZFo) 3 }, (70) 

Recall that in the computation of the PDF tails of the momentum flux 
or zonal flow structure formation, the latter was assumed to be driven by a 
modon in the ITG turbulence. Since this modon is created by the forcing and 
that F(t) = as t — > — oo, we could treat Eqs (66) and (70) as the transition 
amplitudes from an initial state with no fluid motion to a state with different 
values of momentum flux R/Rq or zonal flow formation 4>zf/4>zfo- 

Note that all physical quantities are included in the parameters £i and £ 2 , 
through the ion temperature gradient (rji), density gradient (e n ), temperature 
ratio (r = Tj/T e ), modon size (a), modon speed (U) and wave number (k). 
It is also important to note that £j — > oo (i.e. PDF vanishes) as the the 
forcing disappears (k — > 0); the instanton cannot form and the PDF vanishes 
(P{R) -> and P(0 ZF ) -> 0) without forcing. 
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VI Numerical results of the reduced two field 
model 



We have presented a theory of the PDF tail of structure formation and how 
the PDF tail of momentum flux is modified by the presence of an imposed 
shear flow (Vo). The exponential forms of the two PDF tails are completely 
different, signifying the difference in the physical interpretation. In the case 
of structure formation the PDF tails are found as ~ exp{— ^201^}, while the 
momentum flux PDF tail ~ exp{— £ii? 3//2 }. The origin of these scalings are 
the quadratic non-linearity in the dynamical system (Eqs 1-2). Mathemat- 
ically, the difference in scaling of the PDF tails momentum flux and zonal 
formation comes from the difference in the scaling of the initial condition 
(Fq) with the large parameter A. The spatial structure of the modon is of 
less importance in determining the exponent in the exponential PDF tails 
than the temporal behavior, but changes the overall amplitude through the 
coefficients £1 and £2- Therefore an approximate spatial structure is sufficient 
to determine the exponential scaling whereas the correlation time of / may 
affect the scaling. The spatial structure of the flow is incorporated in the ini- 
tial condition of the flow <fizFo- In this section the parametric dependencies 
of £j will be studied in detail. 

First the PDF tail of momentum flux in ITG turbulence incorporating 
the effects of shear flow is shown in Figure 2. The parameters are r = 2.0, 
e n = 1.0, g t = 1, a = 2, U = 2.0, k = 3000, e = 0.1, k « 1.91 with 
V = 0.0 (Figure 2A), V Q = 5.0 (Figure 2B), V = 10.0 (Figure 2C) and 
V = 15.0 (Figure 2D) for rji = 2.0 (red, dashed line), rji = 4.0 (blue, solid 
line) and 77, = 6.0 (black, dashed-dotted line). When V = 0, the result 
recovers the previous finding in Ref. [Sj. It is clearly shown that the PDF 
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Figure 2A 



Figure 2B 
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Figure 2: The PDF tail of momentum flux in ITG turbulence incorporating 
the effects of shear flow for V = 0.0 (Figure 2A), V = 5.0 (Figure 2B), 
V Q = 10.0 (Figure 2C) and V = 15.0 (Figure 2D) for rj, = 2.0 (red, dashed 
line), r/i = 4.0 (blue, solid line) and r\i = 6.0 (black, dashed-dotted line). 

tail of momentum flux is significantly reduced if a strong shear flow is present 
whereas weak flow can increase the PDF tail. In the equations the flow speed 
and the modon velocity come in as a combination of the form (U — Vo), 
determining the behavior of the resulting PDF tails. When (U—Vq) decreases, 
the PDF tail increases until it eventually decreases. This means that there 
exists a negative value (U — Vo) that gives maximum PDF tail, depending on 
all other parameters. Note that although the Galilean invariance of Eq. (1)- 
(2) may allow us to perform Galilean transformations to find other instanton 
solutions, those solutions correspond to different ground states in the field 
theory and must be discarded since they have nonzero velocity in the infinite 
past [28]. Thus, our instanton solutions are not Galilean invariant |28j . 

Second, the PDF tails of structure formation in ITG turbulence incorpo- 
rating the effects of shear flow is shown in Figure 3. The parameters and the 
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Figure 3: The PDF tail of structure formation in ITG turbulence incorporat- 
ing the effects of shear flow for V = 0.0 (Figure 2A), V = 5.0 (Figure 2B), 
V Q = 10.0 (Figure 2C) and V Q = 15.0 (Figure 2D) for rj, = 2.0 (red, dashed 
line), r/i = 4.0 (blue, solid line) and r\i = 6.0 (black, dashed-dotted line). 

interpretation of the results are the same as those in Figure 2. 



VII Discussion and conclusions 

In summary, this paper presents the first prediction of the PDF tails of zonal 
flow formation from a dynamic system of turbulence and zonal flows. The 
PDF tails of zonal flow is PDF ~ exp{— czf$\f\ (where j — 1, 2 represent 
the different constants depending on the models) which corroborate earlier 
findings [5]- [9]. 

In the first part a general five field model (ni,(/),Ti,T e ,Viu) where a re- 
ductive perturbation method is used to derive dynamical equations for drift 
waves and a zonal flow is studied. Interestingly the predicted PDF tails of 
momentum flux in the five field model of ITG turbulence are Gaussian when 
the feedback of zonal flow on the turbulence is incorporated. This result 
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differs from what has theoretically been found earlier in Ref. [H] and found 
in the reduced two field model of ITG turbulence in Sec. V whereas it is in 
agreement with predictions from non-linear simulations of turbulence [19|- 
[2U] and in gyro-kinetic toroidal simulations [21]- [22] ■ Mathematically it is 
due to the different highest non-linearity between the two models (i.e. cubic 
in the five-field model and quadratic in the two field model). Physically, it is 
because in the two-field model the feedback of zonal flow on the fluctuations 
was not treated self-consistently, which will regulate turbulence, inhibiting 
its growth. 

Note that non-linear models with cubic non-linearities for the turbulence 
fluctuations was obtained perturbatively and thus valid only for small values 
of the amplitude (small 0i). The statistical property of the fluctuations 
differs in systems with cubic and quadratic non-linearities. To elucidate 
these differences it is instructive to discuss a general formula for PDF tails 
of any moment in non-linear systems derived previously [24] . By using the 
instanton method, the PDF of the m-th moment was shown to be 

P(Z) ~ exp{-cZ s }, (72) 

s = 73 

m 

in the case where n is the highest non-linearity in a dynamical system. 

For instance, in the case of the five-field model, n = 3 [Eq. (9)] while in 
the case of the two-field model, n = 2 [Eq. (1)]. For PDFs of momentum 
flux, m = 2 while for PDFs of zonal flow formation m = 1. According to Eqs. 
(72)-(73) there is a significant difference in the PDF tails of the first moment 
in a system with cubic non-linearity where s = 4 and one with a quadratic 
non-linearity where s = 3 [25J. For weak zonal flow Eq. (7) becomes linear 
and we expect Gaussian statistics (s = 2) for the fluctuations. Note that for 
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the five-field model, Eq. (72)- (73) suggest that the PDF tails of zonal flow 
formation are determined by the quadratic non-linearity in the dynamical 
equation for the zonal flow whereas the cubic non-linearity determines the 
PDF tails of momentum flux. 

Finally, in the second part a reduced two-field model (0,7*) that has 
an exact non-linear solution (bipolar vortex soliton) of ITG turbulence and 
zonal flows are studied to compute the PDFs of zonal flow formation and 
momentum flux. One of the important results from the numerical study in 
this reduced two-field model, which is also supported in the five field model, is 
that shear flows can significantly reduce the PDF tails of momentum flux and 
zonal flow formation. Since zonal flows are more likely to be generated further 
from marginal stability, they will then regulate ITG turbulence, leading to 
a self-regulating system. Namely, while ITG turbulence is a state with high 
level of heat flux, it also generates stronger zonal flows that inhibit transport. 
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A Coefficients 

We have balanced the non-linear transport due to the imaginary parts of a 
[^(q;)] (A. 15) with the source term S T i, 

S Tl + 2^^(a)^(0S 5 / 3 )>«-0^0S^) 

- ^wj^ ^r~ 01 ^r^ - (A - 1} 
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Ci = 
x 
+ 
+ 
+ 
+ 

C 2 = 
+ 
+ 

+ 
+ 

c 3 = 
+ 

Cnli = 

X 

+ 

+ 
+ 
+ 
+ 
+ 
+ 



-k s y (l + r)e n (3r 2 en [10e n + 3r(2 - 3^)] 

u{-k y X + ujf + fcjj(25A;Je^[5(l + r)e n 

r(-2 + 3i]i)]\ 2 + 45A; 2 re n [5(l + r)e n 

r(-2 + 3^)]A 2 cu + 27k y r 2 [5(l + r)e„ 

r(-2 + 3^)]AV + r{50e 2 + 15re n [(2 - 3^) + 

9r 2 A[(4-6^) + 5 + 3r)\]}u 3 ), 

fcjL(l + r)e n o;(5A; J/ e n + 3™)(™ 2 {/cJe n [5(-7 

5£>l)e n + 3r(-2 + 3^)] + 30fc y (-l 

D])T6 n u + 9(-l + D 2 ± r 2 to 2 } 

kj[5k 2 y e n [5(l + r)e n + r(-2 + 3^)] 

Qk y r[5(l + r)e n + r(-2 + 3^)]cj 

3r 2 (5 + 3r)w 2 ]), 



^L3P 2 ra + / 2 A; 2 )A; m -4^]), 

A(5e n + 3rA) { ^ W(5 ^ + 3Ta;) 
[k\\(ky{-50£\ - 3r 2 [10e n + 3(2 - 3^)]A 2 

9r 3 (-2 + 3^)A 2 - 10re n [3A 2 + 5e n ( Ve + A)] 

3r 2 [10e n r] e + 3(1 + r)(-2 + 3^)A]cu) 

cu(A; 2 e n [175 en A + 9r 3 (2-3^)A 2 

35re n [3A 2 + 5e n (r] e + A)] + 3r 2 A[3(2 - 3^)A 

5e n (5r) e + 7A)]} + 3/c J/ re n [50e n A 

3r 2 A[(-2 + 10r7 e + ^) + 10A] 

r[10e n (Qr] e + 5A) + 3A(-2 + 3r] t + 10\)]}u 

9r 2 (5e n + rA)[A + r(^ + A)]cu 2 )]}, 
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C N L2 = jz ; — --{k A l 2 T{l +r] i + r\)uj 2 ( y 5k y e n 

(5e n + rA) y 

+ Tu) 2 {10k 2 e n - u[5k y e n (7e n + 3rA) 

+ 3r(5e n + 3rA)cj]}), (A.6) 

Catls = D 2 ± k y l 2 r(l + r)e n o; 3 (5 k y e n + 3ro;) 3 , (A.7) 

C 4 = D*jfeJZr(l + T)e n w 3 (5^e n + 3ro;) 3 , (A.8) 

(1 + r)e n A 2 (5e n + 3rA) 
+ 5re 3 A(677 e - 5D 2 ± X) + 9Z>5r 2 A 4 (l 

+ 2 Vi r\) + L>lre n A 3 [10(l + Vi ) + (22 

- 3r)rA] + e„A 2 {3r[2 + r(2 + 5 Ve - 3r] t ) - 3^] 

+ 5£>i[7(l + 77i) + 2(7 -3r)rA]}), (A.9) 

D 2 = -(3D 2 ± rX 3 (l+r ]l + rX)+5e 3 n [r ]e -(l + r)X] 

+ e 2 A{[7 + r(7 + 10% - 3^) - 3^ + 5L>1(1 + ^)] 

+ r[5Di-3(l + r)]A} 

+ e n A 2 [3r(l + r + r Ve ) + D 2 ± {-7 + 3r)(l + ^ + rA)]) 

x T7 — , (A.10) 
A(l + r)e n ' 1 ; 

^J-k m l 2 (k 2 y D NL1 + D NL2 ) : (A.ll) 

D NL1 = - t T [5e 2 n + (-7 + 3r)e n r-3rX 2 ], (A.12) 
(l + r)e n 

Cmi = 3e„A ! ,»{— }, (A.13) 

*L = (^-A) ^-'"-:!' , (A.14) 



« = ""geA J*" ' ( A - 15 ) 
3r k y 

(A - l6) 

A n = ^{^(1 + ^) + ^}, (A.17) 
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r e n s 
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